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Abstract 



We studied the effect of nearest pattern interaction on a global pattern 
formation in a 2-dimensional space, where patterns are to grow initially from 
a noise in the presence of a periodic supply of energy. Although our approach 
is general, we found that this study is relevant in particular to the pattern 
formation on a periodically vibrated granular layer, as it gives a unified per- 
spective of the experimentally observed pattern dynamics such as oscillon and 
stripe formations, skew-varicose and crossroll instabilities, and also a kink for- 
mation and decoration. 
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Typeset using REVTgX 



Patterns arising in parametrically forced extended media such as Faraday systems, vi- 
brated granular layers, and electrically driven liquid crystals have drawn interests in recent 
years [0] . From a dynamic point of view, it is noteworthy that recent experimental studies of 
the patterns arising in a vertically vibrated thin granular layer have led to the discovery of 
the so-called oscillons, a fluidlike localized nonlinear wave [0,^. Following the discovery of 
oscillons, patterns with different dynamic characteristics in a periodically vibrated granular 
layer were also reported. These are (i) the patterns featuring the skew-varicose and crossroU 
instabilities in a fluid convection and (ii) the kinks; the interfaces or fronts that separate 
domains of different values of phases 

An understanding of the fundamental nonlinear dynamical behavior of granular materials 
remains a serious challenge Several theoretical approaches, including molecular dy- 

namics simulations, order parameter equations, and hydrodynamic-type models, have been 
proposed to describe this phenomenology. Direct molecular dynamics simulations repro- 
duced a majority of patterns observed in experiments and many features of the bifurcation 
diagram, although until now have not yielded oscillons and interfaces. Hydrodynamic and 
phenomenological models @-|n| reproduced certain experimental features; however, neither 



of them offered a systematic description of the whole rich variety of the observed phenom- 



ena 



p!2| . As it is, it is often argued that a general local, continuum description of granular 
media, analogous to hydrodynamics, is unlikely to exist. 

Earlier, we proposed a nearest pattern interaction model to capture the essential phe- 
nomenology of the birth of oscillons. As this model was proposed before the experimental 
reports of the other types of patterns (i) and (ii) mentioned above, it remains primitive. In 
this study we would hke to generalize this model to capture also the dynamic phenomenology 
of the patterns (i) and (ii). 

We start by presenting and briefly outlining the nearest pattern interaction model which 



we proposed earlier to study oscillons. For more details, readers may refer to the Ref. [13 



M : K+^{r) = F[K{r) + aG{AK{m. n = l,2,.... (1) 
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with 

AM^ ^ is / ^r'rih-^'^'^ - hnm'r' (2) 

, where r' runs over entire space but Wp^ is a weighting function being 1 if i? < |r' — r | < 3i? 
but otherwise, and AS* denote the area of the concentric region where W = 1. We set 
hn{r') to be coupled with its neighbors in terms of the introduced averaged fluctuation 



Ahn{r), which effectively allows the interaction only with its nearest patterns (refer to Fig. 
1(a)). This model employs a discrete mapping for time evolution. For a periodically forced 
oscillator, the discrete time step often corresponds to the period of forcing. 

We assume that patterns are to grow from a noise, subject to a statistical constraint 
that a rise(/i > 0) or a fall(/i < 0) from the reference flat surface(/i = 0) is equally likely 
over time. Such constraint satisfies the mass conservation. To implement such assumption 
in the model, we took the following odd function as our functional form of F (See Fig. 1(b)). 

h for \h\ < ho, 

ho for h > ho, (3) 
—Hq for h < —ho. 

With such functional form of F, the noise does not grow and remains as a fixed point if 
their is no nonlinear coupling with neighbors, i.e., if G = 0. 

For the following specific form of G (See Fig. 1(c)), we learned in our earlier study that 
the model is capable of producing the oscillons, and by varying the coupling constant a, the 
model could capture the formation of stripes, oscillon chains and oscillon lattice, etc |]T3 . 



G{Ah) = tan^Ah/{2ho) (4) 

The significance associated with other forms of F and G was a question in our earlier 
study which we could not answer in the absence of other forms of patterns (i) and (ii). With 
more patterns available now, we may proceed to generalize the model as follows. 

Let's start with the function G. The function G is meant to incorporate the interaction 
between neighboring excitations, so we can change the nonlinearity of the interaction by 



using another function for G. Earlier, we had chosen highly nonlinear function Eq.(4) for G, 
to give sufficient hysteresis for stable oscillon structures. What if we choose a less nonlinear 
function for G? In such a case, the model exhibits skew-varicose and crossroU instabilities, 
while the stability region of oscillons becomes narrower. For example, using a simple linear 
function G{Ah) = Ah, and setting a = 1.5, we obtain the pattern dynamics as shown in 
Fig. 2. We have checked many forms of G such as piecewise linear functions and simple 
nonlinear polynomials. In these cases, we obtain in general the patterns as featured in Fig. 
2 as long as G is odd and monotonically increasing. However, these instabilities do not 
occur if G is highly nonlinear as in the case of Eq.(4), for the stability boundary becomes 
too wide. We note that the stability of convection rolls in a fluid heated from below is 
limited by these secondary instabilities. Recently, an experimental study has shown that 
such pattern dynamics are also found in granular systems 0]. 

Next, we turn to the role of the function F. We notice immediately that by changing 
the slope of the function F, we can make the system a bistable system. For that purpose, 
we introduce a new parameter f3 as follows. 

ph for p\h\ < ho, 
Fih) = < ho for f3h > ho, (5) 
—ho for Ph < —ho- 

For j3 > 1, h = becomes a unstable equilibrium, while h = +ho and h = —ho become 
stable equilibria(see Fig. 3.), making the system bistable. Conceivably, as there are now 
two stable attractors separated by an unstable one, the domain may be divided into two 
parts, one with high elevation corresponding to h = +ho and the other with low depletion 
corresponding to h = —ho, separated by a phase front, namely, a kink. The model now 
contains two control parameters, a, (3, which have a definite role respectively. That is, a 
measures the degree of nonlinear coupling with the neighbors, leading into self organization, 
and (3 determines the degree of nonequilibriumness of the system. 

For f3 > 1, and for the tangent function of Eq. (4) for G, the results are summarized in 
Fig. 4. When a is below a certain critical value, ad, indicated by the filled circles in Fig. 



4, the map synchronizes /I's with their neighbors, so flats occur. In this case, however, there 
now emerge two possible phases of flats which are separated by an interface called a "kink" 
0. Figure 5(a) shows a typical kink solution of our model. As we increase a above Oci, 
the interfaces are decorated by periodic undulations as shown in Fig. 5(b,c). We found that 
the areas embedded within the region of the opposite phase as in Fig. 5(b) are unstable 
and shrinks, while fronts which extend to the wall(as in Fig. 5(c)) are stable. We point out 
that the undulated kinks shown in Fig. 5(b), (c) resemble the the decorated kinks reported 
recently in shaken granular beds In our model, decorated parts oscillate at / = 1/2 

while flats are / = 1. But in experiments, both flats and decorations are oscillating at 
/ = 1/2. This is because, in real experiments, flat regions undergo a period doubling while 
a front separating these regions does not 0. 

Associated with the kinks, the present model further shows the followings. First, as a 
exceeds ac2, corresponding to the region denoted by the filled squares in Fig. 4, the decorated 
kinks begin to be elongated to generate stripes. Fig. 6(a) shows the decorations elongating 
into stripes as a is raised to the stripe region(Q; = 16). The same phenomenology is observed 
in experiments [§,0]- Second, the decoration can be converted into oscillon chains when we 



decrease /? below unity, as in Fig. 6(b,c). This point implies that oscillon chains can emerge 
from the decorated kinks when vibrational acceleration is decreased abruptly to the oscillon 
region, and we have confirmed that this phenomenology is indeed observed in a laboratory 
experiment |1^ . This point is noteworthy because the undulation in our model is not merely 



a simple curvature of interfaces, but related with oscillons [|T5[. This aspect is similar to the 
stripes turning into a chain of oscillons as conformed experimentally. 

In conclusion, we have constructed a dynamical model to understand the key mechanism 
for the pattern formation in a parametrically forced spatially extended medium. The nearest 
pattern interaction model here contains two control parameters a and /5; /3 controls the 
number of stable equilibrium states of the system, and a monitors the nearest pattern 
interaction. Despite its extreme mathematical simplicity, our model turns out to provide a 
unified perspective of the various patterns observed experimentally in periodically vibrated 



media. 
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Figure Captions 

Fig. 1. (a) A nearest pattern interaction approximation. The field at P is coupled with the 
fields at other locations only within the shaded region. Here R represents the domain 
size of an excitation. The functions (b) F and (c) G used in Ref. are presented. 

Fig. 2. Secondary instabilities; (a) Skew- varicose instability; (b) CrossroU instability. 

Fig. 3. The function F {(3 > 1). The slope /3, which changes the local dynamics of /i, 
infiuences the overall pattern dynamics drastically. See the text for details. 

Fig. 4. Phase diagram in general. The open circles and squares are to indicate the region 



of hysteresis. Conditions for numerical simulations are the same with Ref. |]13[. See 
the text for details. 

Fig. 5. Various types of kinks; (a) linear kinks (a = 8.0,/? = 1.5); (b,c) Decorated kinks 
(a = 13.0, /3= 1.5). 

Fig. 6. (a) Kinks under kink to stripe transition (a = 16.0, /5 = 1.5). (b,c) Oscillon chains 
from decorated kinks in Fig. 4(b,c) (a = 13.0,/? = 0.95). 
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